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I.  IN  THQDDCTI3 K 


Because  of  threats  from  highly  maneuverable  high 
performance  aircrafts  and  the  need  for  increase  standoff 
ranges,  major  improvements  are  needed  in  guidance  and 
control  capabilities  of  missiles. 

The  high  maneuverability  of  targets,  has  led  to  defense 
missiles  capable  of  develop  higher  lift  accelerations  and  to 
more  complex  control  lavs,  able  to  improve  performance  over 
well  know  lavs  as  proportional  navegation. 

In  order  to  accomplish  these  r.av  requirements  with  large 
standoff  ranges,  propulsion  systems  using  airbraathing 
engines  has  been  studied  and  developed  in  recent  years. 

The  advent  of  airbreathing  engines  leads  natural  to  a 
consideration  of  ban k-to-turn  missiles  in  order  to  minimize 
the  angle  of  attack  of  the  inlets. 

The  necessity  of  more  complex  control  laws,  leads  in  a 
general  way,  to  the  application  of  modern  control  and  esti¬ 
mation  theory,  since  more  complete  informations  cf  the 
states  of  missile  and  target  are  necessary  than  those  states 
informed  by  sensors  commonly  in  use  in  missiles  today.  This 
leads  to  the  use  of  a  airborne  computer. 

The  present  work  adr esses  the  design  and  evaluation  of  a 
optimal  digital  control  for  application  to  terminal  guidance 
in  a  tank- tc- turn  missiles. 

One  continuous  two  dimensional  model  was  adopted,  in  the 
following  form: 

X(t)=  A  (t )  x(t>*  B  (t)  u ( t }  ♦  E  g  (1.1) 


1 1 


■>  v>  kV.v.w.v'.  >  v  -.  'Tv  •  ■ 


where  the  effect  of  gravity  appears  explicitly  in  the  third 
term  on  the  righ  hand  side  of  expression  1.1. 

After  the  development  of  a  equivalent  discrete  model, 
the  optimal  control  problem  has  been  solved, using  a  modified 
Ricatti  equation  due  to  the  existence  of  the  third  term 
representing  to  the  gravity  effect. 

Next,  several  analysis  has  been  made  in  order  to  check 
the  effects  of  small  and  large  roll  excursions,  the  effect 
of  the  sample  rate  on  the  system,  and  the  effect  of  the 
initial  pitch  angle,  in  order  to  check  the  validity  cf  such 
two  dimensional  model,  when  applied  in  some  scenarios  of 
in  te  r  st . 


II.  MODEL  OP  THE  STSTEH 


&.  IHTBODUCTIOH 

In  -he  present  work  the  problem  of  terminal  guidance  for 
long  range,  banfc-to-turn  missiles  with  ramjet  engines ,  using 
a  digitalized  system  has  been  investigated. 

The  model  developed  in  reference  1  is  used  as  the  bass 
for  this  work.  After  the  digitalization  of  that  model,  an 
optimal  control  law  was  developed. 

B.  ASSUMPTIONS 

Keeping  the  same  assumptions  as  in  ref. 1, one  has: 

The  missil  is  limited  to  -2g's  and  +15  g's  of  commanded 
acceleration  in  the  pitch  plane, with  zero  lag. Also  its  yaw 
auto  pilot  has  zero  lag,  yaw  regulator  maintains  zero 
sideslip. 

Missile  thrust  exactly  cancels  drag. 

The  angle  of  attack  is  assumed  to  be  very  small, which 
leads  to  the  commanded  acceleration  acting  normal  tc  the 
velocity  vector. 

The  missile  will  not  have  to  roll  through  a  large 
angle. (Further  considerations  will  give  to  this  at  the  end 
of  the  derivation  of  the  control  law)  . 

C.  THE  COHTIHUOUS  MODEL 

Using  the  same  reference  frames  as  in  ref.  1, one  assumes: 

-Body  frame  with  x„  axis  parallel  tc  the  longitudinal 
missile  axis,  positive  yb  axis  out  of  the  left  wing  ,and 
positive  zb  axis  upward,  (see  fig. 2.1) 
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*  -1*  ■ 


-Fligth  path  axis  with  xf  axis  parallel  to  the  velocity 
vector,  positive  z*  axis  pointing  upwards  and  axis 

pointing  to  the  left,  (see  fig.  2.2) 

In  fig  2.1  and  2.2,  the  angles  $  and  6  are  the  Euler iar. 
roll  and  pitch  angles. 

The  state  vector  is  given  as 

[_  T 

Yf  Y'  Aty  zf  zf  Atz  Aft  J  (2.1) 

where  y„  and  z,  are  the  components  of  the  relative  target 
position,  yp  and  z,  are  the  relative  target  velocity, Atv  and 
Atz  are  the  components  of  target  acceleration , which  is  expo¬ 
nentially  decaying  with  a  time  constant  5  . 


A  ft  -  </>-  /i  (2.2) 

where  is  the  initial  roll  angle  (at  t  =  0)  . 

The  control  vector  is  given  as: 


u 


(2.3) 


where  Ac  is  the  commanded  acceleration  and  Pc  is  the 


commanded  roll  rate. 

The  nonlinear  plant  equation  is 


A  ty  /■  Ac  s  av  <p 
•Aty/S 


(2.5) 


A  ty  f-  Ac  s/a/ 

0 

-Aty/S 

0 

• 

zr 

0 

Atz  -  Ac  cosyi 

-  j  cos  9 

- Atz  /<? 

0 

Pc 

0 

- 

whare  g  is  gravity's  acceleration  and  9  is  the  pitch  angl 
as  seen  in  fig.  2.2 

Linearizing  and  setting 


G  =  e  g 


(2.6) 


one  has 


X  = 


y  f 

Aty+A'c  (cosfrfo)  £0 
- Aty/S 

At  a  +  Ai(cos  6)  &<f 


0 

0 

0 

0 

O 

0 

o 


(2.7) 


5 


h 


2  =  0  0  0  0  “  CCS  &  0 


in  eqn.  2.1 , vs  have  set 


cosfi  -  cos  ( ^  +Afi)  -cos  cosJj^-sin^  sin  Ap 

sin  P  -  sin  { pa  +  Ap)  =sin^5  cos^^cosj^  sin  Ap 

and  expanded  in  Aft  which  is  considered  small. 

Now  assuming  that  t!c ,  which  is  actually  t! 
control  Ac,  can  be  expressed  in  the  fora  of: 


Ac  =  A  co 


[•-  vVl 


with  Aco  =  Ac  at  t  =  0,  one 


x=Ax  +  Bu+Eg 


where 


o  /  o  o 


OO/O 


0  0  ’'/&  O 

OOOO 


OOOO 


OOOOO 


/-  — 

1  77 


Aco  /-  —  sm  d 
T7  to 


(2. 12) 


0  0 

J'O  <f>0  o 

0  0 

5-0  € 

-  cos  4.  0 


o  0 

0  1 


E 


^0  0  3  0  -cos  6 


0 


(2.13) 


g  =  g 


(2.  14) 


D.  THE  DISCHETE  MODEL 

1 •  Introduction 

With  the  introduction  cf  a  digital  computer  to 
control  the  continuous-tine  system,  one  has  to  have  some  kind 
of  interface  in  order  to  take  care  of  the  communication 
between  the  discrete  and  continuous- time  systems. In  this 
case  it  will  be  considered,  A-tc-D  and  D-to-A  converters  as 
samplers  and  zero-order  holders  as  in  reference  2. 

In  such  case,  considering  tae  system: 


i(t)  =  A  (t)  x  ( t)  ♦  B  (t )  u  ( t)  ♦  E(t>  g  (t) 


(2.  15) 
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on 9  can  write  the  state  of  the  system  it  time  t(k  +  1)  a~  : 

*  > -  <*-*,  «v  <*pb<v  °<v  (2.16) 

r  **♦'  $r 

♦/^(t4l  ,*))  E(^)  d*|  g(tk) 

where  ^Z$(t,to)  is  the  transiction  matrix  of  the  system  repre¬ 
sented  by  eqn.  2.13. 

Furthermore,  we  will  considered  that  the  sampling 
instants  are  equally  spaced ,  or: 


:  -  t  =  T 

k<  i  K 


t  =  kT  ♦  r 


so  one  can  replace 
=  kT 

thus , 


(kT+T)  =  0  (kT+1)  x(kT) 


, *ur 

♦  fi  {  kI*T,/>p  3(*p  d*j  a(kD 

/*r*r 

0  (kT*T,/*|  )  S ( /*]  }  d"|  g  (kl) 


(2.  17) 


(2.18) 


(2.19) 


or  in  a  simplified  notation: 


x  ( k  + 1)  *  \d  (k)  x(k)  ♦  Bj  (k)  u  (k) 


(2.23) 


♦  E4OO  g<*) 

wh  2  r  e  , 


A4  (k)  =  fi  (kr  +  T  ,k!) 


r .  -  . 


2 .  Calculation  of  the  Maurices  A  (k)  ,  B(k)  a 


It  is  straightforward  to  show,  using  rha 
of  rhe  matrix  A(t),  that  the  rransiction  marrix 
2-19  is  : 


fi(Jcn-T)  = 


/ 

T 

Aj 

a 

0 

0 

0 

Hr 

0 

f 

Hi 

0 

0 

0 

Ae^I.T 

0 

0 

iT/* 

0 

0 

0 

0 

0 

0 

0 

f 

T 

Hr 

0 

0 

0 

0 

/ 

Hr 

0 

0 

0 

0 

0 

e5* 

0 

(2.21) 


(2.  22) 


(2.23) 


d  E(k) 

sparseness 
of  equation 


(2.  24) 


where: 


AdL.  =Ad  =  5  T  -  5  2  <1-e'%  ) 

For  the  calculation  of  the  others  terms  or.e  may  male 
use  of  the  property  of  the  transictioa  matrix  that: 


JAihJd.  ,  1<v  Ai(t  ti, 

dt. 


so. 


±AA.LkJiTj-kJ2 

d(kT*T) 

For  A  (2,7)  : 

d  Ad„  (kT,T,kT) 


A  (kT+T)  Ad  (kT+T,kn 


d(kT-rT) 

iL^i4/-0i7iL_!sLL  - 

d(kTi-T) 


*  Aat  (kT+T)  =  ACO 


A. 

Ti 


:osK 


ACO  1- 

Ti 


■  ^s/9 


/kT*T  1 

ACO  cos^  /  [i-  ackx^n 

*<*r 

4  **] 


Ad  =  ACO  COS 
A  r 


J  kt-r 


For  A  (1,7)  : 


d  Ad.,(kT*T,  wt) 

- L-  (kT*-T,  kT) 

d(kT*r) 


Ad  (kT+T,  kT) 


fu 


kT.T 


4  (kT«-T) 


•'  -■  '■■.If  F.l'.l.f- 


TS]  »« 


Doing  the  same  process  for  Ad  (5 ,7)  and  Ad(«,7)  one 


has: 


ld..,  ■ ico  sln^.  [  TZ-  11 J 

A  =  Aco  sin/o  [  T  -  ~Yri'~  ]  T*  ] 

For  the  derivation  of  the  marrix  Bd  (k) 
according  to  egn.  2-22 

j6  (kT+T)  B(  j  )  =  Ad  (kT+r,^)  B(ij) 

where. 


AdfkT*Tty  ) 


/ 

Am 

U* 

Am 

0 

0 

0  Am 

1 

0 

/ 

A/n 

'i.i 

0 

0 

0  Am 

■*.r 

(9 

0 

Am 

0 

0 

0  o 

0 

0 

0 

/ 

Aon  4  m  A» 

'*■*  '*.T 

0 

0 

0 

0 

/ 

A\> 

0 

0 

0 

0 

0 

*1"  o 

0 

0 

0 

0 

0 

0  1 

one  needs 


where  A*j  represents  Aj(kT+r,ij  ),  and 
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AdOcT*T,,j  )&(*))  r 


■V!^  Ai«.> 

•“>  A  A"i„ 

0  0 

0  0 


whsre 

kT+T-  1  =  KnUs 


k\r  [ 

kT+T- 

iillll—  .  «  *  .Jjj.  ] 

2  Ti  1  2  Ti  J 

Aco 

cos/i 

T  (kr+I 

’»  -  [“/fr*'  ]  <“♦*> 

r2 

- 

- 

T  «\  ♦ 

COS 

A 

A\r 

j^T(kr*D  -  (kr+T) 

T2 

- 

- 

T  <T|  ♦ 

-Vri— ]  T*l]A=’ 

sin 

A 

n 

<w 

•* 

kT+r- 

Hl.'.lt.  -  *  *  4--  1 

2  Ti  1  z  Ti  J 

Aco 

sinX 

where 


/kT  *T 

(kT+l-^)  d*j  sin ^  Sin^ 


\,u  '  T  C0SA 

/  kT* T 

3d«*  I  M,.,1! 

*<T 

which  after  some  algebric  work  has  been  found 


B  (k) 


f  _Tl_  _  _  (kT*T_/_kT_  _  (kt/_  J 

[  7  JTi  ^  r;  6Ti  J 


ACC  cos 


,fcT*  T 


BdU(k)  =  /  %d1 


which  can  be  found  to  be: 
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■  l  7"  *  V-  T<J 


ACO  COS 


'*• ,.*»■  [  7l 
[  7- 


•  -r  ["I?,’”]  T*J 400  Sir-A 

CjcT  +  t)*  (  kT  +  T )  *  .  (kT)  1  .  / 

- - --  > - —  kT  - - -  A  co  siap 

3  Ti  2  t;  <Ti  I 


In  the  same  way: 


0(*T.T,y)  £(*))  = 


-  T  COS  9 

-  cos  © 


thus  E  (k)  is  equal  to 

3  =  £  0  0  0  0  - j-  cos©  -r  cos©  0  0  j 

where  we  have  considered  0  as  a  constant  angle.  (Further 
comments  cn  this  after  development  of  the  control  law). 

Notice  that  througout  this  work  ,the  commanded 
acceleration  has  been  considered  an  unknown  and  the  assump¬ 
tion  has  been  made  that  it  will  be  in  the  form  of  eqn.2.9. 
One  will  need  in  further  developments  to  consider  the 
control  Ac  as  a  known  Ac  (k)  ,  which  will  baa  constant  between 
kT  and  kT+T.  With  such  assumptions,  the  discrete  represen¬ 
tation  of  the  system  is  easily  found  to  be  : 


x  (k  +  1 )  *  A  (k)  x  (k)  *  E(k)  u  ()c)  ♦  15  g 

where 


/ 

r 

4</ 

i.i 

0 
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Ad  =  3(1-  e'T/*)  -  Ad 

2$  J 

T*  ,  / 

K.r=  ~Y  C0S^ 


Ad47=  T  AC  C0s/o 
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=  -r-  Ac  sm/a 


Ad 


«./ 


Ad 


T  Ac  sin 


(2.  26) 


(2.27) 
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III.  THE  OPTIMAL  COST HOLLER 


* 


i.  DERIVATION  OF  THE  OPTIMAL  CONTROLLER 

In  order  to  have  a  suitable  guidance  law  ;o  implement 
the  control  commands,  we  will  minimizing  the  following 
performance  index: 


where  x  (N)  is  the  final  state  at  t=Ti 

As  we  want  to  minimize  the  final  miss  distance, the 
weighting  matrix  W(N)  is  taken  as 


and  Q(k)  is  a  two  by  two  positive  definite  symmetric 
weighting  matrix  to  be  chosen. 

In  the  derivation  of  the  solution,  reference  3  has  been 
followed  keeping  in  mind  that  the  state  equation  has  th3 
form : 
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x  (k  +  1)  =  a  (k)  x  (k)  ♦  B  (k)  u(k)  ♦  S  j 


or 


\ 


x  (k  +  1)  =  f  (x  (k)  ,u  (k)  ,g) 


whare  the  third  t  e  rm, whir  h  represents  the  s 
gravity  has  baer.  considered  constant. 

Considering  that  the  parfcrmanca  index  is  : 


J  =  /£x(N)J+\  L  (k)  £  x(k)  ,u  (k)  ,g(k)J 


k  »  O 


we  need  to  find  a  seguence  of  u(k)  that  minimi: 
Adjoin  tha  system  equation  to  J  with  a  mul 


J=0[x(N)]^^  L(k)  £  x  (k»  ,u(k)  ,gj  + 

k.*o 

♦  X  (k+1)  |  fK  £x  (k)  ,u  (k)  ,gj  -  x(k+1)j 

and  defining  a  scalar  sequence  H  (k) 

H(k)=  L(k)  ^  x(k)  ,u  (k)  ,g  J  + 

+  X  (k+1)  £k  £  x  (k)  ,u  ( k)  ,  g  J 

k=0 ,1,2,.... 


f  f  ec 

p.  th 

as  J 
ipli 


n- 1 


one  has: 
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with  boundary  condition 


X(N) 


.  ijL 

c)  x(w) 


(3.12) 


we  obtain  the  minimization  of  the  performance  index. 
In  the  present  case  wa  have: 


J=  y  xT(N)  W(N)  x  (N)  »  y  uT(k)  2  ( Jc)  u(k) 
♦  AT(k+1)  [  A  (k)  x  <k)+B(k)  u(k) +E  3-x(k  +  1)  J 


(3.  13) 


and  H  (k)  : 


H(k)  =  y  uT(k)  2  (k)  u  (k)  > 


(3.  14) 


XT  (k*1) 


A  (k)  x  (k)  ♦  3{k)  a  (ki  ♦  E  g 


then  in  order  to  minimize  J: 


"**  s  aT(k)  Q  (k)  ♦  XT(kf1)  B(k)  =3 


(3. 15) 
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or  considering  that  Q(k)=Q(k) 
Q  <k)  u  (k)  =  -  BT  (k)  >  (k+  1) 


XT(m  = 

^  xfk) 


>!(k)  =  >T(k  +1)  A  (k) 


(3.16) 


£<N)  =  xT(N)  M(N) 


(3. 17) 


Notice  rhat  we  are  not  weighting  the  states  in  the 
performance  index,  except  the  last  state. A  more  general  form 
could  be  obtained,  with  all  states  being  wei gthing  ,  if  we 
change  the  eqn  3.16  to: 


£(k)  =  >>  (k  +1)  A  (k)  ♦  xT(k)  Ml  (k) 


(3.18) 


where  Ml  (k)  is  the  weighting  matrix  of  the  states. 

with  equations  3.15,3.16  and  3.17  one  is  able  to  find 
the  sequence  of  u(k)  that  will  give  the  minima  controls. 

Such  set  of  equatics  can  be  solved  by  the  sweep  method 
as  in  ref. 2  . 

Me  will  look  for  a  solution  of  the  form: 


u  (k)  =-F  (k)  x  (k)  -FG  (k)g(k) 


(3.19) 
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what  aaans  that  the  commanded  acceleration  and  rc 
will  have  a  correction  due  to  the  effect  of  gravity, 
Placing: 

A(k)  a  S(k)  x(k)  ♦  SG(k)  g  (k) 


from  eqn.  3.15 

Q  (k)  u  (k)  «-BT  (k)  Js<k+1 )  x(k+1)  +3  3  (k+1)  g  (k  +  1)] 
from  egn.  3.3 


3(k)  u(k)  =-BT(k)  S  (k-el )  £  A  (k)  x(k)«- 

♦  B  (k)  u(k)  ♦  3  g]  -  BT(k)  SG(k  +  1)  g(k  +  1) 
so 

[q  (k)  ♦  BT  (k)  S  (k  +  1)  B(k)J  u  <k)  = 

-BT  (k)  S  (k+1)  A  (k)  x  (k)  -  BT  (k)  5  (k  +  1 )  E  g  (k) - 
-BT  (k)  SG(k+1)  g  (k  +  1) 

considering  that  g  is  a  constant 

=  -  Tq  (k)  +  BT  (k)  S  (k+1)  B  (k)l 


.22) 


u(k) 


(k)  S  (k+  1 )  A  (k)  X(K)  + 

♦  |  BT  (k)  S  (k  +  1)  E  +  Br  (k)  Sg(k+1)J  gj 

so : 

u  (k)  =-F (k)  x(k)  -FG  (k)  g 

where 

F  =  [q  (k)  ♦  B T  (k)  S  (k+ 1)  B(k)J  .  (3.2 

•  BT(k)  S  (k+1)  A  (k) 


FG 


*  £q  (k)  +  BT  (k)  S  (k+1)  B  (k) 


(3. 


[bt  (k)  S  (k  +  1)  E  ♦  BTSg  (k+1)J 
from  eqn.3.16  and  3.19 

X(k)=AT(k)  A(k+1)=AT(k)  [  S(k+1)  x(k+1)  +SG  (k  +  1)  g  (k+1)J  r 


=AT  (k)  S  (k  ♦  1 )  A  (k)  x  (k)  ♦  A'  (k)  S  (k+1 )  3(k)  u(k) 


+  AT  (k)  S  (k ♦  1 )  Eg  +  AT  (k)  SG(k+1)  g(k+1) 


from  eqn.  3.19 


X(k)  =  AT(k)  S  (k  +  1)  A  (k )  x  (k)  +  AT  (k)  S(k  +  1)  B(k). 


[-F(k)  x(k)  -  FG (k)  gj  +  AT(k)  S(k+1)  Eg  ♦ 
♦  A  SG  (k)  g (k+1) 


so,  as  g  is  a  constant: 

S  (k)  x  (k)  +SG  (k)  g=  1  At  (k)  S  (k+1)  A  (k)  -  A  (k)  S  (k+1  )B  (k)  F(k)1 


+\k  (k)  S  (k+ 1}  E  -  A  (k+1)  8(k)  FG (k) 
+  A  (k)  SG  (k+1)  ]  g 


S(k)=A  (k)  S  (k+ 1)  A  (k)  -A  (k)  S  (k+ 1)  B  (k)  F  (k) 


(3.  25) 


SG  (k)  =AT  (k)  S  (k  +  1)  2-AT  (k)  S  (k+1)  B  (k|  FS(k) 


(3.  26) 


♦  A  (k)  SG(k+1) 


Thase  equations,  3-25r  3.26,  3.23  and  3.24  can  be  solved 
backwards  with  the  final  condition: 

S  (N)*W(N) 

SG  (N)  =0 

Nctics  that  this  satisfies  our  previous  boundary  condi¬ 
tion  in  eqn.  3.17  where: 

(E)«i(R)Z  (N) 

S  (N)  X  (N)  =W  (N)  X  (N) 


S(N)=»(N) 

B.  EXTEHTIO*  OF  THE  HO DEL  FOE  LASGE  BOLL  ANGLES 

Op  to  this  point,  one  has  to  taka  into  account  that 
throughout  the  development  of  this  work, the  angle  has  been 
considered  small,  it  is  necessary  to  relax  this  restriction. 
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In  order  to  do  that,  the  system  has  been  broker,  in  two 
blocks  as  in  figure  3.1. 

The  first  block  is  a  representation  of  the  algorithm 
which  will  calculate  the  optimal  commands.  The  algorithm  has 
contain  with  itself  an  exact  model  of  the  system  or  missile. 
The  model  of  the  missile  is  initialized  from  the  information 
on  the  initial  states,  the  initial  input  command  Aco  and 
initial  roll  angle  (at  t=0  ) ;  and  computes  the  optimal  gains 
and  further  the  optimal  commands  which  will  be  feed  to  the 
missile. 

The  method  adopted  in  computing  the  optimal  commands  is 
more  easily  understood  if  one  considers  figure  3.2. 

In  figure  3.2,  the  lines  numbered  as  0  in  the  graph  for 
Ac  and  for  Pc  are  the  optimal  commands  for  a  given  initial 
roll  angle(^0).  Lines  number  1  are  the  commands  for  a  second 
initial  roll  angle  (£$,)  larger  than^,  and  so  on.  Thus  in 
figure  3.2  one  has  a  family  of  optimal  commands  for  any 
initial  roll  angle. 

Notice  that  the  upper  line  of  the  graph  of  Ac  represents 
the  accelerations  of  a  missile  which  had  at  t=0  a  correct 
initial  roll  angle  in  order  to  hit  the  target  with  no 
comands  in  roll  rate. 

In  the  present  method  the  computer  performs  the  calcula¬ 
tion  of  the  commands  only  for  the  firs-  step  of  time  and 
then  feeds  these  commands  to  the  missile.  The  missile  is 
then  driven  to  the  next  state  (x(k+1|)  and  feed-backs  to  the 
computer  the  information  on  the  roil  angle  at  that  step.  The 
roll  angle  feed-backs  from  the  missile  is  considered  by  the 
algoritm  as  the  initial  roll  angle  at  t=0  and  the  next 
commands  are  calculated. Not  ice  that  at  this  second  step  the 
algorithm  will  feed  to  the  system  the  second  command  (at 
t=t,  )  .  This  process  is  them  repeated  until  t  is  equal  to 
the  intercepted  time. 


It  is  important  to  realize  that  with  this  method  of 
calculation,  since  the  algorithm  was  developed  with  "-.he 
assumption  of  small  roll  excursions  soma  error  is  expected 
due  to  the  fact  that  in  computing  the  gains  by  solving  a 
Ricatti  equation  backwards,  as  has  been  done, it  is  necessary 
to  update  the  system  from  t=0  to  t=ri  at  each  step,  and  in 
this  process  the  roll  angle  is  not  small. Notice  however  that 
we  are  applying  the  commands  only  in  one  step, and  if  one 
expects  that  the  roll  rate  will  decrease  to  zero, as  we  are 
increasing  in  time,  the  variation  of  the  roll  angle  will 
tend  to  decrease, sc , we  can  expected  that  the  error  will 
decrease  as  the  time  increases. 

Another  important  point  tc  be  studied  is  how  the  missile 
itself  (second  block  in  fig. 3. 2  )  has  to  be  implemented  in 
order  to  be  valid  for  small  and  large  roll  excursions. 

It  is  considered  that  one  has  the  perfect  knowledge  of 
the  commands,  thus  the  missile  is  modeled  as  a  state  vari¬ 
able  system  as  in  eqn  2.26  with  the  initial  roll  angle  being 
update  at  each  step.  In  this  way  the  system  will  take  the 
initial  roll  angle  as  the  summation  of  all  previous  initial 
roll  angles,  (see  fig  3.4) 

Prom  the  original  variables  one  has  for  the  state  6<f>  , 
the  following: 

A^k+1)  =  A^(k)  +  Pc  (k)  T  (3.27) 


which  is  shew  in  f igg.  3.3,  where  the  initial  roll  angle  is 
kept  constant  and  is  update  each  str-p. 

Notice  however  that  considering  large  roll  excursions 
(see  fig.  3.4)  and  keeping  in  mini  that  the  angle  ft  has 
been  defined  as 
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tha  expression  3.27  is  not  valid,  since  in  modeling  the 
system  it  was  assumed  that  would  ba  samll. 

This  leads  to  a  change  in  the  expression  for  the  varia¬ 
tion  of  as  in  fig.  3.4.  In  fig.  3.4,  the  angle  is 

update  at  each  step,  so  ona  has: 

Afifr*  11=  Pc  (k)  T  (3.  28) 


By  consideration  of  the  equation  2.27  it  can  be  seen 
that  by  setting  tha  element  A {7, 7)  to  zero,  one  can  obtaim 
equation  3.28. 

Thus  the  missile  model  for  large  roll  excursions  will  be 
represented  by  tha  following  equation. 


x(k  +  1)  =  A  (k)  x  (k)  ♦  B(k)  u  (k)  ♦  3  g 
where 
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where  A(7,7}=0  ,  for  large  roll  excursions 
A  (7,7)*i  ,for  small  roll  excursions 


Y  s,n& 

T  Sif\  to 


T2 

■  _  C  03  j 

2 

T  cos  0 

'o 


*T*S  y 

—  Ac  cos  0 
6 

T*  , 

—  Ac  cos)£ 


T3  / 

-j*  Ac  s'*#, 

-I-*  Ac  t0 


(3.  30) 


-  -If  cos  0 

z 

.  r  cos  0 


(3.31) 


We  will  redefine  the  states  x,  ,  as  the  relative  posi¬ 
tion,  x2  as  the  relative  velocity  in  the  Y  direction  and  x3 
as  the  target  acceleration  in  Y  direction.  The  states  xv  , 
xs  and  x6  has  the  same  meaning,  out  in  the  Z  direction  and 
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xr  is  .  In  -this  representation  of  the  system,  x,  (k)  is 
equal  to  the  component  of  the  miss  distance  along  the  Y 
direction  and: 


x,  (fc+i)=  y„(k+i)=  y,  (k)  +  y,  (*)  t  + 


(3.32) 


♦  [c  I-  52  (1-  e  ^  )J  Aty  (k)  Ac(k)  cos^  4^(k) 

Tl  ,  T*  j 

♦  -j~  smyo  Ac(k)  +  -j-  Ac(ki  c^spo  Pc(k) 

The  first  three  terms  in  the  SH3  of  eqn.3.32,  are  easily 
seen  as  the  contri buit ion  to  the  miss  distance  of  respec¬ 
tively  the  previous  miss  distance,  relative  velocity  and 
target  acceleration.  The  following  two  terms  represents  the 
contribution  of  the  commanded  acceleration  and  the  last 
term  represents  the  effect  of  coupled  Ac  and  Pc,  and  tends 
to  be  small  due  to  the  cube  of  the  sample  period. 

For  the  component  of  miss  distance  in  the  Z  direction, 
one  has: 

(k+1)=  Zp  (k+  1)  =  z  .  (k)  +  zF  (k)  T* 

♦  [s  T-  S  2  (1-  e*%  )  Atz  (k)  ♦  -j-  Ac  (k)  sin^4^(k)  - 

-  -j-  cos^Ac(k)+  Ac(ki  sinjtf  Pc  (k)  -  -j-  ccs  0  g 


where  its  ‘eris  have  the  same  physical  meaning  as  in  the 
expression  for  x  (k+1),  with  the  effect  of  the  gravity  added 
to  the  expression. 

Since  one  can  notice  that  in  tae  representation  of  the 
miss  distance  in  Y  direction  appear  two  terms  as  a  function 
cf  cos  ,  and  in  the  rep resentaction  of  the  miss  distance 
in  Z  direction  appear  two  terms  as  function  of  sin  jZ^  ,  it  is 


interesting  to  verify  that  the  fourth  term  in  the  5HS  o 
both  expressions  acts  like  a  correction  for  the  firry  term 
Befering  to  fig.  3.4b,  one  can  see  that  at  any  step  of  rime 
the  commanded  acceleration  is  actually, 

Ac  ccs^-^Ac 

and  considering  small  angles: 

4  Ac=  Ac  cos  j6a  -  Ac  cos(^e  +  &$)  * 

=  Ac  cos £  -  Ac  £cos^Z$  4^-  sin ^4 sin  AjtJ  = 

=  Ac  sin)4 

The  same  idea  can  be  applied  to  the  expression  fo 
x  (k  + 1)  . 

The  terms  x^  and  xf ,  represent  the  relative  velocity 
an d  are: 

x2  (k+1)=  y  (k)  ♦  s  (1"  -  ^  )  Aty(k)  *■  Ac(k)  cos^  T 

jl 

♦  Ac (k)  s in/4  T  +  ~J~  Ac(k)  cosjz£  Pc(k) 

xs  (k+1)=  z„(k)  ♦  5  (1-  e'^6  )  Atz(k)  *  Ac(k)  sin j6o  T  - 

-  Ac<k)  ccs  T  +  ~2~  Ac(k)  sin^  Pc(k)  -  T  cos  0  g 

Where  the  two  first  terms  in  the  RHS  represents  the  effec 
of  the  velocity  and  acceleration  at  a  previous  step,  and  th 
ether  three  terms  has  the  same  meaning  as  previously  stated 
The  terms  x,  and  x4  are  the  target  accelera tions,  i 
this  model  being  exponentially  decaying. 

1 .  Effects  on  the  Miss  Distant 2  of  the  Sxtsnticn  of  th 
Model 

In  previous  subsection,  a  attention  of  the  model  fo 
large  roll  excurssions  has  been  performed.  Notice  t^t  ther 
are  two  models  of  the  system  being  used.  The  first  one,  use 


in  the  algorithm  is  valid  only  for  samll  roll  excursions, 
and  a  second  model,  valid  for  small  and  large  roll  excur¬ 
sions  used  as  a  representation  of  the  missile. 

The  algorithm  with  the  first  model,  as  explained 
before,  is  initialized  at  each  step  with  the  actual  roll 
angle  of  the  missil,  and  performs  the  calculation  of  the 
comands. 

In  order  to  check  the  effect  of  the  extention  of  the 
model  on  the  miss  distance,  one  can  define  a  ideal  initial 
roll  angle  (  jrfo,atAk  ),  as  the  roll  angle  at  t=0  in  order  to 
have  the  comanded  acceleration  vector  pointing  to  the 
projected  final  target*  s  position.  This  means  '•hat  the 
missile  would  not  have  to  roll  to  hit  the  target  (see 
fig3.5  ),  thus  the  commanded  roll  rate  calculated  by  the 
algorithm  will  be  equal  to  zero.  This  implies  that  from  that 
point  ahead,  the  roll  angle  is  constant  and  equal 
to  r  and  that  the  time  history  of  the  control  Ac  will 

be  a  straight  line. 

The  fact  that  the  roll  angle  will  tend  to  this  limit 
deserves  an  investigation.  Notice  that,  as  show  in  fig.  3.2, 
at  the  moment  the  missile  reachs  its  maximum  roll  angle,  the 
control  Ac  will  be  the  required  acceleration  to  hit  the 
target  if  the  missile  initial  roll  angle  was  •  This 

means  that  the  control  Ac  computed  would  be  correct  only  if 
the  missile  would  have  turned  immediatly  to  this  angle. 

For  very  small  angles,  the  previous  comment  would  be 
acceptable,  but  in  a  normal  situation,  as  showed  in 
figure3.5,  the  missile  will  only  reach  the  ideal  initial 
roll  angle  after  some  time,  and  due  to  the  vertical  compo¬ 
nent  of  the  acceleration,  when  this  occurs  the  missile  would 
be  in  a  position  above  the  ideal  trajectory,  which  means 
that  it  would  follow  a  coarse  parallel  to  the  ideal  trajec¬ 
tory  to  intercept. 
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Ore  can  see  that  such  problem  will  lead  to  a  large 
miss  distance,  in  the  case  that  the  target  acceleration  it 
is  not  small.  Thus,  some  correction  is  necessary  in  order  to 
improve  the  missile  performance. 

Figure  3.5b  shows  the  missil  at  some  point  of  its 
trajectory  where  it  has  reach  its  maximum  roll  angle,  thus 
it  is  at  a  parallel  course  with  its  ideal  trajectory  to 
intercept.  kt  this  point,  since  all  the  states  are  known, 
it  is  possible  to  recompute  a  new  ideal  .  So,  if  at  this 
point  the  computer  is  feed  with  with  the  states  at  this 
point  it  will  compute  the  commands  in  order  to  drive  the 
missil  to  the  new  0^d<al  ,  which  will  introduce  the  desired 
correction.  Notice  that  such  correction  can  be  made  during 
all  the  flight,  from  t=0,  to  t=Ti. 

In  the  present  work  this  method  has  been  accom¬ 
plished  by  feeding-back  to  the  computer  the  roll  angle  and 
all  the  states  of  the  missile.  with  this  information  the 
computer  is  able  to  perform  the  calculation  of  the  corrected 
commands  at  each  step  of  time.  howeveras  the  states  are 
being  updated,  it  is  also  necessary  to  update  the  time  to 
intercept,  which  has  been  done  using  the  time  to  gc  to 
intercept,  or: 

Ti  (k)  =  Ti-  k  T 

where  Ti  is  the  nominal  time  to  intercept. 

C.  SIHtJLATION 

In  order  to  keep  the  same  assumptions  as  reference  1, 
the  matrix  Q  (weighting  the  control)  has  been  put  as 
suggested  in  such  reference  or: 
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7V‘.v, 


r 


0 


(3.33) 


b, 


with 

b,  =5.78  103 

b^=5.0  meters 

Five  different  cases  were  ran: 

Case  0,  tested  with  one  simple  model  valid  only  for 
small  roll  angles  begins  with  missil  and  target  on  parallel 
courses  to  the  inertial  x  axis,  the  target  100  meters  above 
the  missile  and  with  an  evasive  manoevar  exponentially 
decaying  with  time  constant  of  23  seconds.  The  initial 
acceleration  of  target  was  -.5  g's  in  y  direction,  (.’ee  fig. 
3.6) 

Case  1,  with  the  same  scenarios  as  case  0,  but  was  run 
using  the  algorithm  for  large  roll  excursions. 

Case  2, is  the  same  scenario  as  in  case  1,  except  that 
the  initial  acceleration  of  target  was  -1.0  g's  in  -.he  y 
direction.  (Same  as  case  1  in  ref.  1i  . 

Case  3,  the  same  as  case2,with  target  acceleration  of  -4 
g's.  (  Same  as  case  2  in  the  reference  1) 

Case  4,same  as  case  3  but  with  target  at  initial  posi¬ 
tion  600  meters  bellow  the  missile.  (Same  as  case  3  in  refer¬ 
ence  1.) 

D.  COMMENTS  AND  CONCLUSIONS 

1 .  Results 

In  case  0r  the  missile  bagins  its  trajectory 
commanding  26.5  m/sec2  and  the  time  history  of  the  control 
Ac  follows  exactly  a  straight  line  in  the  form  suggested  in 
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ref.  1,  as  seen  in  fig. 3.  7The  control  Pc,  baginir.g  a*  .35 
rad  per  second,  is  decaying  and  ceachs  zero  at  k=80  (see 
fig.  3.8).  Figures  3.9  and  3.10  show  the  miss  distance,  where 
one  can  see  that  the  missile  is  crossing  the  target  with  a 
CG-to-CG  distance  of  1.5  meters. 

Figure3. 11,  shows  the  time  history  for  the  roll 
angle,  which  as  expect,  reachs  a  constant  value,  with  the 
missile  crossing  target  at  t=Ti,  with  a  bank  of  .44  radians. 

In  case  1, which  was  run  with  the  model  for  large 
roll  excursions,  the  missile  kept  the  same  Aco,  but  there  is 
a  very  small  increase  in  further  commanded  acceleration  in 
order  to  correct  the  effect  of  the  roll  angle  on  the 
vertical  component  of  the  control  Ac,  as  seen  in  fig.  3.  12, 
and  table  I. 

Figure  3.13  shows  that  tie  roll  rate  decreases 
almost  as  before,  and  the  final  roll  angle  is  .42  rds,  as 
shown  in  fig.3.15  and  table  I.  Sefering  to  figure  3.14  there 
is  a  change  in  the  final  miss  distance,  which  is  better  than 
case  0,  due  to  a  improvement  in  its  Z  component  (see  fable 
I)  .  This  results  in  a  final  CG-to-C3  miss  distance  of  .65 
meters. 

In  case  2,  the  missile  has  the  same  Ac  at  t=0,  but 
with  the  correction  for  roll  angle  being  increased  due  to 
the  increase  of  target  acceleration  (see  fig. 3. 16),  the 
commanded  Ac  reachs  a  larger  peak  value- 

The  initial  value  of  the  commanded  roll  rate  is  .69 
radians,  which  is  larger  than  case  1,  due  to  the  increase  in 
the  target  acceleration.  In  figure  3.18,  there  is  no 
noticeable  change  in  the  shape  of  the  curve  for  Z  direction, 
and  in  the  Y  direction  the  final  distance  is  about  the  same 
as  in  case  1.  Tne  CG-to-C3  distance  at  t=Ti  is  .73  meters. 
The  larger  rcLl  rate  leads  to  larger  roll  angles  as  seen  in 
fig.  3.19,  where  the  final  roll  angLa  is  .72  radians. 
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The  effect  of  target  acceleration  can  be  easily  seer, 
in  case  3,  where  one  can  see  that  with  the  sane  Aco,  the 
accelerations  are  largelly  increased  from  this  point,  and 

the  missil  begins  its  trajectory  with  very  high  roll 

rate  (see  figures  3.20  and3.  21)  .  There  is  a  change  in  the  Z 
component  of  the  miss  distance,  that  decreases  its  final 
value  to  .05  meters,  but  now  the  miss  distance  in  Y  direc¬ 
tion  is  made  worse  as  shown  in  figure  3.22,  which  leads  to  a 
final  CS-to-C3  distance  of  1.5  meters.  In  fig. 3. 23  one  can 
see  that  the  missile  crosses  the  target  with  a  bank  angle  of 
1.26  radians. 

In  case  4,  due  to  the  position  of  target  600  meters 
under  the  missile,  the  initial  commanded  acceleration  is 

negative  and  reachs  the  limit  of  -2  g's.  The  initial  roll 

rata  begins  at  a  smaller  value  than  in  case  3  but  increases 
during  the  initial  part  of  the  flight  reaching  its  peak 
value  at  1.75  seconds  whan  again  as  in  the  previous  cases 
begins  to  decay.  As  the  missile  banks  tc  roll  angles  larger 
than  90  degrees,  the  acceleration  goes  to  positive  values, 
as  seen  in  figures  3.24  and  3.25  Figure  3.26,  shows  the 
worse  case  among  these  in  respect  to  the  miss  distance, 
mainly  in  the  Z  direction,  and  in  the  final  cg-to-cg 
distance,  which  is  equal  to  4.43  meters.  Also  the  final 
roll  angle  of  3.0  radians  is  the  largest  among  all  these 
cases,  as  seen  in  fig.  3.27. 

2 .  Comments 

Defining  the  projected  zero  effort  miss  distance 
(Z2H )  as  the  miss  distance  the  missile  cross  the  target  with 
no  commands.  It  can  be  calculated  at  t=0  as  the  initial 
distance  between  target  and  missile  plus  the  miss  distance 
due  to  the  gravity  or: 
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ZEM  =  z  (0)  ♦  -j~  g  t* 


(3.34) 


which  is  equal  in  all  the  three  first  cases. 

In  the  first  three  cases,  the  initial  missile’s 
commanded  acceleration  is  the  same  as  seen  in  table  3.1. 
Considering  that  the  control  Ac  necessary  to  correct  the 
initial  miss  distance  can  be  calculated  as: 


Ac 


t<tu. 


i 

l 


26.7 


(3.  35) 


which  is  close  to  the  initial  control  Ac. 

This  suggests  that  the  initial  Ac,  would  be  that  one 
necessary  to  correct  the  initial  ZEM  in  Z  direction,  which 
agrees  with  reference  1.  Notice  however  that  in  all  cases 
the  initial  Ac  is  less  than  the  calculated  value  of  Aco, 
which  agrees  with  the  previous  statement  that  some  error  was 
expected  in  the  initial  part  of  the  computations. 

Considering  the  ideal  initial  roll  angle  as  defined 
before,  one  has: 


(3.36) 

From  table  I,  and  figures  showing  roil  angles,  in 
can  be  verified  that  the  missile  is  banking  to  reach  angles 
larger  than  •  in  order  to  correct  its  trajectory  to 

hit  the  target. 

Therefore,  in  all  cases,  the  missile  begins  its 
trajectory  with  an  Aco  (discussed  before)  ,  and  a  roll  rate 


deal 


=  tan*1 


z£  My 

ze 
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which  is  proportional  to  the  target  accelaration  in  ^he  Y 
direction.  As  the  missile  rolls  at  decreasing  roll  rates, 
the  commanded  acceleration  is  changed  in  order  for  to 
compensate  the  effect  of  the  roll  angle  on  its  Z  direction 
component. At  some  time  when  the  control  Pc  is  zero  or  near 
zero,  the  control  Ac  begins  to  follow  a  linear  law,  as 
suggested  in  eqn.  2.7: 


Ac=  Aco 


[’-  tV] 


Notice  however  that  the  term  Aco  in  this  equation  is  no  more 
the  actual  initial  commanded  acceleration,  but  that  one  the 
missile  would  have  if  its  initial  roll  angle  was  equal  to 
the  final 

Such  behavior  defines  a  boundary  in  the  control  Ac 
which  is  clearly  seen  in  3.29,  where  the  commanded  accelera¬ 
tion  is  bounded  by  the  curve  of  the  control  Ac  the  missile 
would  have  if  its  initial  fa  was  equal  to  ^®ldeBl(°r»  if  no 
commanded  roll  was  necessary  to  reach  the  target)  . 

Althoug  the  missile  commands  roll  angles  larger  than 
the  ideal  initial  rcll  angle,  it  is  possible  to  do  a 
prediction-  with  no  no  computer  work-  of  an  aprox  ima  tior.  for 
the  maximum  acceleration  the  missile  would  experience  during 
its  fligth,  as  folio  wing  (see  fig.  3.32): 

i  (3.37) 

cos 

aiW«ol 

where  fa  ,  cames  from  eqn.  3.50,  and  Ac  is  a  streiaht  line 
as  in  fig  3.32. 

Other  interesting  point  is  than  the  missile  is 
commanding  no  reach  roll  angles  about  twice  of  fa  i#(etl  when 
the  target  is  at  small  accelerations,  and  when  at  large 
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accelerations,  the  missile  is  making  a  small  correction  on 
its  roll  angle,  as  shown  in  table  I. 

One  can  see  from  the  figures  showing  miss  distance, 
that  the  relative  position  of  missile  to  target  is  about  the 
same  in  all  three  initial  cases.  Thus,  the  new  $l<ieai  at 
each  point  is  the  same.  Is  the  0oldea\  computed  at  t=0  is 
smaller  when  the  target  is  at  small  accelerations,  the 
correction  has  to  be  larger  in  order  to  hit  ohe  target,  as 
seen  in  figure  3.30. 

In  table  I,  one  can  see  the  final  miss  distances, 
the  final  miss  distance  eg  to  eg  and  the  final  roll  angle. 
Such  results  show  that  with  this  digitalized  model,  good 
results  has  been  obtained. 


Ac(0) 


Figure  3.2  Variation  of  Commands  with  Initial  Boll  Angle. 
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Figure  3.3  Variation  of  Soil  Angle  -  Snail  Angles 


Figure  3.U  Variation  of 


•:v 


.t; 


M  »l.l  •\mSmlm.  \  t.’.  -  &  UjkA  3  -  3 


Figure  3,5b 


CASE  1  target  acceleratfon  -.5  g 


200  4 


100  4  target  tntttal  poettton-  CASE  1 


-loo  4 


-200  4 


eteefle  tntttal  poettton 


200  4 


100  4  target  tntttal  poettton-  CASE  2  and  3 


case  2 

target  acceleration  -1  g 


-too  4 


-200  4 


-300  4 


-400  4 


-500  4 


etsetle  tntttal  posttfon-CASE  2,3  and  4 

case  3  and  4 

target  acceleration  -4  g 


-000  4  target  tntttal  poettton  -  CASE  4 


Figure  3.6  Scenarios  foe  Simulation. 
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ACCELERAT  1 0N-riETER5/SQ .  SEC 

5.0  7.5  10.0  12.5  15.0  17.5  20.0  22.5  25.0  27.5  30.0 


CA5E  0 

INITIAL  TARGET  ACCELERATION-. 5  G 
INITIAL  TARGET  POSITION-  100  M 
SAMPLE  PERIOD-0.05  SEC 


Figure  3.7  Coaaanded  Acceleration-  Casa  0. 
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RATE-  RDS/SEC 


CASE  0 

INITIAL  TARGET  ACCELERATION— .5  G 
INITIAL  TARGET  POSITION-  IOO  « 
SAMPLE  PERI0D-Q.Q5  SEC 


INITIAL  TARGET  ACCELERATION- >.5 
INITIAL  TARGET  POSITION-  100  M 
SAMPLE  PER I 00-0. 05  SEC 


CRSE  0 

INITIAL  TARGET  ACCELERATION— .5  G 
INITIAL  TARGET  POSITION-  100  M 
SAMPLE  PER I 00-0. 05  SEC 


ACCELERATION  -  HETERS/SQ.SEC 


1ST  CASE 

INITIAL  TARGET  ACCELERATION-  -.5  G 
INITIAL  TARGET  POSIT ION- 100  M 
SAMPLE 'PER 100-0. OS  SEC 


62 


1ST  CRSE 

INITIAL  TARGET  ACCELERATION-  -.5  G 
INITIAL  TARGET  POSITION-lOO  M 
SAMPLE: PER 10D-0. 05  SEC 


Figure  3.13  Coananded  RoLl  Rate-  Case  1. 
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MISS  DISTANCE-  METERS 

-100.0  -ao.o  -60.0  -40.0  *20.0  0.0  20.0  40.0  60.0  60.0 
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ANGLE-  RDS 


ACCELERATION  -  METERS/SQ.SEC 

0.0  2.5  5.0  7.5  10.0  12.5  15.0  17.5  20.0  22.5  25.0  27.5  30.0 


2ND  CASE 

INI  URL  TARGET  ACCELERATION-  -1. 
INITIAL  TRRGET  POSITION-lOO  M 
SAMPLE  PER  1 0D-0. 05  SEC 


2ND  CRSE 

INITIAL  TARGET  ACCELERATION-  -}.  G 
INITIAL  TARGET  POSITION-IOQ  M 

sample: per i od-o. os  sec 


Figure  3.19  Holl  Angle-  Case  2 


ACCELERATION  -  METERS/SQ.5EC 


RATE-  RDS/SEC 
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SflMPLC  PtRiaD-O.CS  SEX' 


30.0  48.0  50.0  60.0 

:  :  if  : 


ROLL  ANGLE-  RDS 

0.25  0.50  0.75  1.00  1.25  1.50  1.75  2.00  2.25  2.50  2.75  3.00 


3RD  CASE 

INITIAL  TARGET  ACCELERATION-  -4.  G 
INITIAL  TARGET  PQSITION-IOO  14 
SAMPLE  PERI0D-0.Q5  SEC 


4TH  CASE 

INITIAL  TARGET  ACCELERATION-  -4.  G 
INITIAL  TARGET  POSITION— 600  tt 
SAMPLE  PER 100-0. 05  SEC 


4TH  CASE 

INITIAL  TARGET  ACCELERATION-  -4.  G 
INITIAL  TARGET  POSITION— BOO  M 
SAMPLE  PER I OD-O. 05  SEC 


ACUU 


Figure  3.28  corrected  Bodel. 
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FTg.  3.30  b  Small  target  accelerations 


Figure  3-30  Corrections  on  the  Roll  ingle. 
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17.  ASALJSTS  OF  Gil NS,  SAMPLE  BATE  AND  PITCH  ANGLE 


A.  ANALISYS  OF  THE  GUNS 

In  chapter  3  has  been  developed  a  solution  for  the 
optimal  control  of  a  systam  as: 

x (k+1 ) =  A  (k)  x  (k)  +  E  (k)  u  (k)  +  E  g 

It  would  be  interesting  to  chack  if  tha  optimal  gains 
reach  steady  stats, but  at  the  moment  that  the  extension  for 
large  roll  excursions  has  been  introduced,  and  the  system  is 
being  feed  with  optimal  commands  which  ara  varying  each  step 
of  time,  such  idaa  can  not  be  applyed.  However  one  can  do 
such  check  in  the  modal  for  small  roll  excursions,  which  is 
acctually  used  to  compute  the  optimal  commands. 

Doing  this,  one  has  the  time  history  of  those  gains  as 
in  4.1  and  4.2 

Ona  can  notice  in  fig.  4,1  that  tha  gain  FS(2,1)  associ¬ 
ated  with  the  effect  of  gravity  has  no  effect  on  the 
ccmanded  roll  rate,  and  that  FG(1,1),  as  shown  in  fig.  4.2, 
has  a  large  affect  on  the  commanded  acceleration. 
Furthermore,  this  gain  reachs  stsady  stats  very  fast,  thus 
it  can  be  assumed  that  the  gain  FG  will  be  egual  tc  the 
steady  stata  value  during  all  time. 

From  eqn.  3. 19,  and  assuming  steady  state: 


u  ( k)  =  -F  ( k)  x  (k)  -  FG  (k)  g 


(h.i) 


and  substituting  g: 


u  (k)  =  -F  (k)  x (k)  -  C 


(4.2) 


where  the  second  term  in  the  right  hand  side  is  a  constant, 
and  its  value  is  exactly  equal  to  the  value  of  the  commanded 
acceleration  necessary  at  t  =0  to  correct  the  gravity  fall  of 
the  missil  (or  to  correct  the  initial  ZEM  due  to  gravity)  . 

It  might  be  supposed  that  one  could  solve  the  optimal 
control  problem  for  the  system  represented  by  eqn.  4.),  just 
considering  one  reduced  system  represented  by: 


x  (k+1) =  A  (k)  x  (k)  ♦  B (k)  u  (k) 


with  a  bias  in  the  control  as: 


u  (k)  =  -F  (k)  x  (k)  *■  C 


(4.3) 


(4.4) 


But  as  shewed  in  the  following  analysis,  this  is  not 
possible. 

The  constant  term  in  the  right  land  side  of  egn.  4.4  is 
calculated  as  follows: 
from  egn.: 

initial  ZEM  due  to  gravity=  ZEM  =  -  -  g  t 2 


_ _ 

[•r  '  vs  ]  cos/° 


and  the  gains  F(k)  are  calculated  using  a  Ricatti  equation 
as  usual. 


Case  5  was  tested  using  the  above  considerations. 


and 


using  the  same  scenario  of  case  3. 

Figures4.3  and4.4,  show  the  time  history  of  the 
controls.  One  can  see  that  the  commanded  acceleration  has 
begun  an  same  values  as  in  case  3,  bun  the  commanded  accel¬ 
eration  reachs  a  peak  considerable  higher,  them  dcreases 
does  not  following  a  linear  law,  with  a  final  at  14.7  meters 
per  second  sguare,  being  this  terminal  value  due  to  the 
constant  term  representing  the  effect  of  the  gravity. 

Refering  to  fig.  4.4,  the  commanded  roll  rate  begins  at 
a  same  value  as  in  case  3,  but  as  the  control  Ac  is  too 
high,  it  reachs  negative  values,  going  to  zero  almost  at  the 
end  of  the  running  time.  This  behaviour  of  the  control  leads 
to  a  large  miss  distance  as  seen  in  figure  4.5,  and  tablell. 

Figure  4.6  shows  the  time  history  of  the  roll  angle, 
which  rises  to  values  close  to  1.5  radians.  As  the  accelera¬ 
tion  at  this  point  is  larger  than  the  correct  value,  the  the 
corrections  are  excessive  and  the  roll  angledecreases  at  the 
end  of  the  running  time  to  the  value  of  .12  radians. 

Thus,  one  can  see  that  the  gain  due  to  the  gravity's 
acceleration  can  net  be  replaced  by  its  steady  state 
value. This  kind  of  simplification  can  thus  not  be  done  in 
the  system  being  studied. 

B.  EFFECT  OF  THE  S  AMPLE  RATE 

Throughout  all  the  simulatios  a  sample  period  of  .05 
seconds  has  been  used. In  this  section  a  brief  study  of  the 
effect  of  the  change  of  this  sample  rate  is  given. 

Two-,  best  c.ases  have  been  selected  to  ilustrate  the 
effect  of  the  sample  period. 

The  first  case,  case  6,  has  been  run  with  a  sample 
period  of  .1  seconds  and  corsits  of  the  same  scenario  as 
case  3. 
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As  one  can  see  in  figar@4.7  and  taolelll,  there  is  no 
noticable  change  in  the  commanded  acceleration,  but  the 
commanded  roll  rate  begins  at  smaller  value  than  in  case  3, 
as  is  seen  in  fig. 4. 8  This  initial  decrease  in  roll  rare, 
leads  to  a  large  miss  distance  in  ?  direction  as  shown  in 
fig. 4.9,  and  to  a  small  value  of  roll  angle  (see  fig. 4.  10). 

The  second  case,  case  7, has  a  smple  period  of  0.025 
seconds.  There  is  nc  noticeable  change  in  the  time  history 
of  the  control  Ac  as  shown  in  fig. 4.  11.  The  commanded  roll 
rate  begins  at  a  higher  value  than  in  case  3  as  shows  fig. 
4.12,  which  leads  to  a  final  miss  distance  in  I  direction  of 
-2.22  meters  and  in  Z  direction  smaller  than  case  6  (see 
fig4.13).  Figure4.14  shows  that  the  final  roll  angle  is 
increased  and  the  missile  cross  the  target  with  1.28  radians 
and  with  a  CG-to-CG  distance  of  2.5  meters. 

In  both  cases,  the  miss  distance  was  increased  over  the 
nominal  value  obtained  with  a  sampling  rate  of  0.05  seconds. 
Thus,  it  would  appear  that  there  is  an  optimal  value  for  the 
sampling  rate,  which  may  be  conected  with  the  geometry  of 
the  scenario  and  with  time  to  go. 

C.  EFFECT  OF  THE  INITIAL  PITCH  ANGLE 

It  is  important  at  this  srep  no  remenber  that  througout 
this  work  has  we  have  been  discussing  a  dimensional  model, 
where  there  is  no  information  on  the  X  coordinate,  so  it  is 
impossible  rc  verify  the  behavior  of  the  pitch  angle. 

In  rhis  work,  since  in  all  the  previous  scenarios  the 
initial  angle  0  was  equal  to  zero,  this  value  has  been  kept 
as  a  constant  during  all  time,  and  considering  that  whirhout 
any  information  of  a  third  dimension  it  was  not  possible  to 
correct  the  time  to  intercept,  this  time  was  also  kept 
constant  and  equal  tc  the  nominal  value  of  5  seconds. 


Notice  that  this  assumption  is  likely  to  be  correct  if 
one  has  the  horizontal  initial  distance  from  target  to 
missile  compared  with  the  initial  vertical  distance  between 
target  and  missile  large  enough  in  order  to  have 
small  angles. 

The  question  that  rises  is,  how  could  this  pitch  angle 
affect  the  system  if  it  was  net  small? 

As  seen  in  fig. 4. 15,  the  missile  velocity  in  the  X 
direction  would  be: 

~  cos8  ♦  g  cos©  sin©  t  (4.5) 

which  has  an  effect  not  only  from  the  pitch  angle,  repre¬ 
sented  by  the  cos  9  ,  but  also  from  the  gravity's  accelera¬ 
tion,  which  will  leads  to  a  different  time  to  intercept. 

Considering  the  same  physical  scenario  as  in  case  4,  but 
changing  the  initial  pitch  angle,  in  order  to  have  the 
missile  pointing  to  the  target  (see  fig. 4. 16),  and  keeping 
the  missiles  velocity  of  10  00  m/sec  in  the  X  direction,  one 
has  a  ccmpletly  different  geometry  of  the  problem  as  seen 
from  the  flight  path  reference  frame. 

With  this  new  situation  (see  fig.  4.  16)  ,  case  8  has  been 
run.  Figures  4.17  and  4.18  show  that  now  the  missile  is 
comanding  large  positive  accelerations,  and  the  roll  rate  at 
the  begining  of  the  flight  is  too  high,  going  to  zero  in  a 
very  small  period  of  time.  The  miss  distance  as  seen  in 
fig. 4.19  are  increased  in  the  initial  parr  of  the  fligth  and 
as  the  missile  corrects  its  trajectory  it  is  decreased  to 
reach  a  final  CG-to-CG  miss  distance  of  2.17  m.  The  roll 
angle,  due  to  the  large  control  Pc  is  oscilatory  in  the 
begining  and  becames  constant  with  a  value  of  .57  radians 
(see  fig. 4. 20  and  table  IV)  . 
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Notice  that  the  high  values  of  acceleration  needed  are 
in  some  part  due  to  a  vertical  component  of  target‘s 
velocity,  which  is  seen  from  the  flight  path  frame  as  the 
target  was  manouvering  in  the  Z  direction  with  constant 
velocity.  These  large  accelerations  leads  to  roll  rates  too 
large  for  the  physical  integrity  of  the  missile.  This  means 
that  although  althoug  the  good  results  obtained,  if  compared 
with  case  4,  they  are  not  practical. 

In  order  to  get  rid  of  the  vertical  manouever  of  the 
target,  case  9  has  been  run.  In  this  case,  the  scenario  is 
the  same  as  before  with  the  target  also  pointing  down,  with 
the  same  pitch  angle  as  the  missile,  and  has  a  X  velocity 
equal  to  the  previous  case  (sea  fig. 4. 16). 

Prom  fig. 4. 21,  one  can  see  that  the  decrease  in  the 
control  Ac  is  substantial  if  compared  with  case  8,  but  the 
commanded  roll  rate  is  still  too  large  as  shows  fig.  4.22  The 
time  history  of  the  miss  distance  is  batter,  resulting  in  a 
final  CG-to-CG  distance  of  about  50*  of  case  8  (see 
fig.  4.23).  The  roll  angle  is  not  oscilatory  as  seen  in 
fig. 4.24  and  the  missile  cross  the  target  with  1.3  radians 
in  roll,  (see  tablalV). 

The  results  obtained  in  the  two  previous  cases,  suggests 
that  the  algorithm  developed  in  this  work  could  be  readly 
applied  to  air-to-surf ace  missiles.  In  the  latter,  the 
scenario  would  be  favorable  to  the  missile  than  in  either 
previous  cases,  since  one  can  consider  that  the  target  could 
be  essencially  stationary  in  comparison  with  the  missile 
speed . 

Case  10  has  been  run  with  this  assumption,  and  the  scen¬ 
ario  as  in  fig.  4.16.  The  target  is  with  zero  velocity  and 
acceleration,  and  the  missile  begins  its  flight  600  meters 
above  it,  with  the  same  initial  pitch  angle  as  before. 

Fig.  4. 25  shows  the  time  history  of  the  commanded  accel¬ 
eration,  where  one  can  see  that  as  there  is  no  roll  rate  to 
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be  commanded,  the  acceleration  is  following  a  straight  line 
with  very  resoanable  values.  The  miss  distance  is  shew  in 
figure4.27,  which  shows  the  final  CG-to-CG  distance  of  .31 
meters. 

Based  in  the  results  of  these  tests,  one  can  see  that 
there  will  be  some  effect  of  the  pitch  angle  on  the  miss 
distance,  not  only  due  to  its  effect  on  the  time  to  inter¬ 
cept,  but  also  because  at  the  moment  that  there  is  a  pitch 
angle  different  from  zero,  even  if  the  target  is  keeping  its 
flight  level,  in  the  flight  path  frame  a  component  of  the 
target's  velocity  will  show  up  leading  the  missile  to 
command  large  accelerations  and  roll  rates.  Althoug  this 
results  harm  the  performance  in  an  air-to-air  missile,  in 
the  case  of  air -to-surf ace  missiles,  when  the  target  has 
been  considered  with  no  motion,  good  results  has  been 
obtained. 

D.  EFFECT  OF  TIME  TO  INTERCEPT 

In  the  simulation  of  case  1  and  2  in  chapter  3, it  has 
been  observed  that  when  the  target  was  at  small  accelera¬ 
tions,  the  missile  did  larger  corrections  on  its  roll  angle, 
with  respect  to  the  ideal  initial  roll  angle,  than  when  the 
target  was  with  large  accelerations.  One  can  think  that  must 
be  seme  kind  of  compromise  between  the  velocity  rate  of 
target  and  missile  (which  will  reflected  on  the  time  to 
intercept),  and  the  relative  position  between  them,  which 
will  affect  the  miss  distance. 

In  order  to  do  a  brief  analisys  on  this,  case  11  and  12 
has  been  run. 

In  case  11,  the  scenario  of  case  2  has  been  kept,  with 
the  exception  that  the  missile's  velocity  was  change  to  2000 
m/sec.,  which  means  that  li  was  changed  to  2.5  seconds. 
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Figure  4.29  shows  that  the  acceleration  is 
increased  due  to  the  small  time  required  to  correct  the  ZEM, 
and  the  commanded  roll  rate  is  almost  twice  of  case  2  (see 
fig.  4.30).  The  final  miss  distance  is  more  than  four  times 
the  value  obtained  in  case  2,  as  seen  in  fig.  4.31  and  table 
V.  The  final  roll  angle  is  about  half  as  in  case  2,  since 
the  projected  final  position  of  the  targes  in  the  I  direc¬ 
tion  is  less  than  in  case  2  (see  fig.  4.32). 

Case  12  was  run  with  a  scenario  less  favorable  to  the 
missile,  where  all  the  conditions  of  the  previous  case  was 
kept,  except  the  target  position  that  has  been  increased  to 
200  meters  above  the  missile. 

Now,  the  commanded  acceleration  are  much  larger,  with  a 
initial  Aco  of  103  m/sec2,  being  almost  impossible  to  see 
the  difference  of  the  time  history  of  the  acceleration  from 
one  straight  line,  as  shown  in  fig.  4.33  The  commanded  roll 
rate  is  small,  about  the  same  as  in  case  2  (see  fig.  4.34). 
The  change  in  the  miss  distance  is  noticeable,  with  a  final 
CG-to-CG  distance  of  4.8  meters  as  in  table  7  ,  and  figure 
4.35.  The  final  roll  angle  explain  the  shape  of  the  acceler¬ 
ation  curve,  since  with  the  small  roll  angle  as  shown  in 
fig.  4.36,  the  system  is  behaving  as  for  small  roll  excur¬ 
sions. 

Notice  that  from  this  analysis,  one  has  to  realize  that 
there  is  some  kind  of  envelope  where  the  2-D  system  is 
valid.  And  in  order  tc  determine  this  envelope,  one  has  to 
take  in  account  not  cnly  the  geometry  of  the  scenario,  but 
also  the  time  to  intercept,  which  is  determined  not  only  by 
the  relation  of  velocities  of  missile  and  target,  but  the 
pitch  angle  too. 
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Figure  4.6  Roll  Angle-Case  5 
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Figure  4.25  coaaanded  Acceleration-Case  10. 
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Figure  4.31  Hiss  Distance-Case  11 
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The  scope  of  the  present  work  was  the  development  of  an 
optimal  digital  control  to  be  applied  on  a  bank-to  turn 
missile. 

A  two  dimensional  model,  as  sugested  in  reference  1,  was 
adopted.  After  the  digitalization  of  the  continuous  model  it 
was  necessary  to  solve  a  modified  Ricatti  equation  since  in 
the  state  equation  there  was  a  third  term  representing  the 
gravity’s  effect.  The  approach  that  has  been  adopted  is  new, 
and  although  good  results  were  obtained  for  the  scenarios 
considered  in  this  work,  is  necessary  that  the  algorithm  be 
further  tested  and  evaluated  in  similar  problems  due  to  its 
novelty. 

The  optimal  wa3  solved  with  an  initial  restriction  to 
small  angles.  This  condition  was  later  relaxed  so  that  large 
roll  angles  could  be  analyzed. 

It  is  dificult  to  compare  the  present  work  with  previous 
results  since  Stallard  has  indicated  a  mistake  in  his  orig¬ 
inal  paper,  and  further  works  in  this  area  was  net  found. 

However  some  comparison  with  Stallard  work  is  possible. 
The  commanded  acceleration  cf  the  missile  are  such  as  to 
correct  the  ZEM  at  each  point,  this  agrees  with  that  refer¬ 
ence.  There  is  a  proportional  relationship  between  the 
commanded  roll  rate  and  the  commanded  acceleration  ,  and  the 
commanded  rcll  rate  is  proportional  to  the  defined  0 ,dea/  at 
each  point,  which  again  agrees  with  reference  1. 

The  algorithm  developed  in  this  work  requires  extensive 
computation  at  each  step,  and  is  clear  that  some  software 
optimization  will  be  needed. The  motivation  for  considering 
the  constant  steady-state  gain  due  to  gravity  was  to 
decrease  this  computational  burden.  This  approach  however 
resulted  in  unacceptable  miss  distance. 
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Another  paint  of  investigation  than  could  reflect  or.  the 
period  available  to  the  computer  to  perform  its  calculations 
was  a  change  in  the  sample  rat9.  Iwo  different  sample  rates 
were  investigated,  both  lead  to  larger  errors  than  the 
nominal  period  of  .05  seconds.  A  detailed  study  cn  this 
issue  is  'left  as  sugestion  for  future  works,  since  seme 
optimal  value  of  the  sampling  rate  is  clearly  indicated. 

It  is  important  to  keep  in  mind  that  the  model  adopted 
is  two  dimentinal,  while  the  actual  problem  is  three  dimen¬ 
sional,  thus  some  brief  studies  were  conducted  in  order  to 
check  the  region  of  validity  of  the  2-D. 

In  the  analysis  of  the  pitch  angle,  one  car.  see  that  is 
necessary  to  have  small  variations  in  pitch  in  order  to 
approximate  it  as  a  constant.  However,  at  the  moment  than 
this  angle  is  different  from  zero,  as  explained  in  chapter 
4,  it  is  possible  to  have  in  the  flight  path  reference  frame 
a  target  manouver  in  the  Z  direction  that  will  lead  to  large 
acceleration  commands,  leading  the  missile  to  large  miss 
distances,  when  considering  a  movable  target.  when  the 
present  system  was  tested  tested  against  fixed  targets,  the 
results  were  quite  gcod,  this  suggests  the  application  of 
the  model  in  air -to-surf ace  missiles. 

Further  investigation  were  made  on  the  effect  of  time- 
to-go.  As  expected,  decreased  time  to  go,  results  in 
increased  miss  distance.  A  detailed  analysis  of  more 
complex  scenarios  is  needed  in  order  to  properly  define  the 
effect  of  time  to  go. 

Also,  it  would  be  interesting  to  extend  the  model  to 
three  dimensions  and  include  the  effects  of  lags  on  the 
system  in  future  works.  Finally,  in  appendix  A,  the 
computer  model  used  in  this  work  is  enclosed.  Some  improve¬ 
ments  in  this  program  can  be  done,  mainly  in  the  data  intro¬ 
duction,  and  in  some  optimization  of  the  running  time. 
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APPENDIX  A 
FORTH  AN  PROGRAH 


These  appendix  provides  a  listing  of  the  computer 
program  used  in  the  present  study. 

Since  the  routines  used  are  non-IMSL,  and  a  small  change 
to  double-precision  was  necessary,  they  are  also  being 
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